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By means of first-principles calculations within the density functional theory, we study the structural and 
optical properties of codoped ZnO nanowires and compare them with those of the bulk and film. It is disclosed 
that the low negative charged ground states of nitrogen related defects play a key role in the optical absorp- 
tion spectrum tail that narrows the band-gap and enhances the photoelectrochemical response significantly. A 
strategy of uncompensated N, P and Ga codoping in ZnO nanowires is proposed to produce a red-shift of the 
optical absorption spectra further than the exclusive N doping and to get a proper formation energy with a high 
defect concentration and a suppressed recombination rate. In this way, the absorption of the visible light can be 
improved and the photocurrent can be raised. These observations are consistent with the existing experiments, 
which will be helpful to improve the photoelectrochemical responses for the wide-band-gap semiconductors 
especially in water splitting applications. 
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I. Introduction 

ZnO is a semiconductor with wide band-gap, and is stable 
in aqueous solution. One of the most remarkable applications 
of ZnO in photoelectrochemical (PEC) field is that, it can be 
used to split water by solar light to produce hydrogen as clean 
fuel. Besides, ZnO is also a promising candidate of photoan- 
ode material in PEC cells'. For most of wide-gap semicon- 
ductors such as ZnO and Ti02, the band-gap is far beyond the 
threshold of 1 .23 eV for the water splitting, even though the 
overpotential looses require a larger threshold energy about 
1.7 eV^. On the other hand, the wide band-gap semiconduc- 
tors can merely absorb the ultraviolet portion of the solar spec- 
tra. To enhance the efficiency of utilizing the solar energy, a 
key issue is how to narrow the band-gap of semiconductors 
below 2.0 eV^ so that the majority of the visible light can be 
absorbed. To solve this issue, a suitable band edge positioning 
with respect to the potentials of reduction and oxidation reac- 
tions, together with the suppression of charge recombination 
centers, should be considered^. 

There has been a number of studies concerning the above 
considerations. Among them, the donor-acceptor codoping 
in ZnO and Ti02 materials is remarkable. Recent theoretical 
calculations show that, it is possible to narrow the band-gap 
of ZnO bulk by codoping (Ga-nN) or (AIh-N)'*; the band-gaps 
of Ti02 bulk can be narrowed to be 2.26 eV, 1.6 eV and 0.9 
eV for charge compensated codoping (MoH-C)-, noncompen- 
sated codoping (CrH-N) and {N+Cf, respectively. Experimen- 
tally, doping Cu^ and codoping (Ga-nCu)^ in p-type, codop- 
ing (GaH-N)"^ or (A1h-N)'" in n-type ZnO films have been re- 
alized. In their optical absorption spectra, there is a long tail 
from the absorption peak in the ultraviolet range that persists 
in the central visible range, thereby narrowing the band-gap. 
The key factors to enhance the PEC responses by codoping 
Ga and N in ZnO films are believed to include either produc- 
ing an excellent crystallinity, or to narrow the band-gap, or to 
generate a wide depletion layer to suppress the recombination 
ratd2. It is also observed that N doped n-type ZnO nanowires 



(ZnONWs)ii can enhance PEC responses that are similar to 
those in codoped ZnO films. A theoretical study on the O- 
vacancy in ZnONWs suggests that the neutral and charged 
states of defects have crucial contributions to optical absorp- 
tion in the visible range'^. The charged states of O-vacancy in 
Ti02 bulks have also been emphasized'^ '"*. 

In this paper, we shall explore the structural and optical 
properties of uncompensated N, P and Ga codoped ZnONWs 
by means of the first-principles calculations within the density 
functional theory (DFT)i^. We have found that the low nega- 
tive charged states of N-related defects play a key role in pro- 
ducing the long tail of the optical absorption spectrum, which 
narrow the band-gap and enhance the PEC response signifi- 
cantly. It is also observed that the uncompensated N, P and 
Ga codoping can generate a red-shift of the optical absorp- 
tion spectra larger than the exclusive N doping, and a proper 
formation energy with a high defect concentration and a sup- 
pressed recombination rate can be acquired. These observa- 
tions will be useful for the application of ZnONWs in the PEC 
field (e.g. the water splitting). 

This paper is organized as follows. In Sec. II, we shall 
briefly introduce the calculational approach, including the 
generalized gradient approximation plus U (GGAh-U) correc- 
tion, and the extrapolation method to correct the band-gaps. 
In Sec. Ill, the negative charged states and optical absorption 
spectra of N-doped ZnONWs will be explored. In Sec. IV, 
the defect bands and uncompensated N, P and Ga codoping in 
ZnONWs to enhance the PEC response will be discussed. An 
extrapolation method is applied to discuss the relevant exper- 
imental observations. Finally, a summary is given. 

II. Calculational Approach 

The calculations are performed using the DFT based VASP 
code'*, which employs the pseudopotentials of the projec- 
tor augmented wave (PAW)'^, with the GGA in the Perdew- 
Burke-Ernzerhof (PEE) formJ** for exchange correlations. 
The cutoff energy for the plane-wave basis set is 400 eV. The 



2 



convergence criteria for DFT calculations are 10~^ eV of total 
energy and 10^"^ eV/A of force. All structural relaxations and 
formation energies are calculated with spin polarization, with 
the Morkhorst-Packk-points meshes of 4x4x4, 5x5x1 and 
1 X 1 X 7 for ZnO bulk, ZnO film and ZnONW, respectively. 
After a full relaxation, the density of states (DOS) and the 
optical absorption spectra of ZnONW are calculated without 
spin polarization, and the k-points mesh is raised to 1 x 1 x 21. 

To obtain the absorption spectra, the frequency dependent 
dielectric matrices are calculated by VASP after structural re- 
laxation. The real and imaginary parts of the dielectric tensor 
are obtained by the usual Kramers-Kronig transformation and 
a summation over empty states, respectively''', from which 
the energy dependent dielectric function along the axial direc- 
tion of the nanowires and thus the imaginary part of refractive 
index can be obtained. The optical absorption intensity is pro- 
portional to the imaginary part of refractive index. 

It is known that the conventional DFT method always un- 
derestimates the band-gap of semiconductors, and the posi- 
tion of the charged states of the defects with respect to the 
valence band edge also tends to be underestimated. The char- 
acter of the defect bands is a mixture of the valence and con- 
duction bands, which can be reflected by the change of the 
defect bands from GGA to GGAh-U calculations^". Note that 
in the following. Figs. 2-6 are calculated by GGA without 
U correction, while for Fig. 7, the GGAh-U method^' imple- 
mented in VASP code is employed to calculate the optical ab- 
sorption spectra for several key cases. Dudarev's approach^^ 
is adopted, where only the difference of the Coulomb and ex- 
change parameters {U — J) is meaningful. 

Inspired by the extrapolation scheme in previous works 
about ZnO^"-^^, based on the results from both GGA and 
GGAh-U with a small value of {U — J) for d electrons, we 
extrapolate the calculated band-gap of the ground state of N 
doped ultrathin ZnONW to the band-gap of the experimen- 
tally synthesized N doped ZnONW", and then estimate the 
band-gaps and the optical absorption peaks if other defects 
are doped in the same ZnONWs in experiments. The extrapo- 
lation equation that is similar to the one given in Ref. H^is as 
follows: 



E,{D)=E^^^+^{D) + 

{No) ^ EfGA^No) 



E\ 



GGA+U 



[i?ff(iVo) 



jpGGA+U 

a 



(1) 



where Nq represents substitutionally doping one N atom in- 
stead of O atom, and D represents other defects. Any band- 
gap Eg corresponds to a peak in the optical absorption spec- 
trum, no matter how it is obtained either from GGA (Ef'=''^) 
and GGAh-U {Ef^'^+^X or from the experiment iE'g''P). 

Besides the extrapolation method to correct the band-gap, 
it is instructive to note that by incorporating U not only for d 
electrons but also for s or p electrons, the effective corrections 
to the band-gaps and spectra have already been obtained in 
several casesi2ii2i2i. To aim at giving an approximate but intu- 
itive view on the absorption spectra, we shall consider several 
situations for a comparison in the following [as will be seen in 



Figs. 7(a) and (b)], namely, U is incorporated for the d elec- 
trons of Zn, p electrons of O'^, as well as the corresponding 
electrons of their substitutional dopants. 

In a supercell, the formation energy E^ of a charged defect 
D with a net charge q can be calculated through the following 
equation 
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E^iD") - E'°\D'i) - E'°'{pure) - J:,An,fi, + q^l,, (2) 

where and E*°*{pure) are the total energy with and 

without the defects, respectively, An^ is the change of atomic 
number, fii is the reservoir chemical potentials, i — Zn, O, 
N, P, Ga, and /ig is the Fermi level (FL) with respect to the 
valence band edge varying between zero and the band-gap Eg. 

For charged defects, a jellium background is used. How- 
ever, we note that the periodic boundary conditions may 
induce spurious electrostatic interactions between charged 
unit cells, possibly leaving errors in energy from the first- 
principles calculations^^. For an aperiodic three-dimensional 
condensed matter system, Makov and Payne evaluated the to- 
tal energy by^^ 
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where Eo, q, a, L, e, and Q is the desired electrostatic en- 
ergy of the point defect in the dilute limit, the ionic charge 
of the defect, the Madelung constant of the lattice, the lin- 
ear dimension of the supercell, the dielectric constant, and the 
quadrupole moment independent of L, respectively. The L^^ 
term is the Madelung energy for the lattice of point charges 
in the jellium background, and for quasi one-dimensional sys- 
tem such as nanowires, L corresponds to the periodicity of 
the nanowire^**. The L^^ term exists with a quadrupole mo- 
ment and is inversely proportional to the supercell volume^-. 
If the vacuum separation between the nanowires is fixed and 
only the axial length I of supercell changes, the supercell vol- 
ume is proportional to I. More methods to correct the forma- 
tion energy have been demonstrated^*^'^^ for several cases of 
nanowires. 

In the present work, to evaluate the effect of the energy er- 
rors on the structural stability, we model the energy correction 
term similar to the form of Madelung energy. The total energy 
is expressed as 



E'°' = Eo + j, 



(4) 



where C is the parameter and is zero for neutral charged 
nanowires, and I is the axial length of the supercell. Then, 
Eq. (2) is corrected to be 



E^{D''J) = E^{Di) + 



C 



(5) 



E^{D'')=EQ{D\l)-Eo{pure,l)-j:,An,n, + qfi,. (6) 

Eq (-D') is the desired energy of the defect in the dilute limit. 
For a given charged defect D"^ with given /i^ and /ig, Eq{D'^) 
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TABLE I: The formation energies of neutral charged ZnONW:N 
under O- and N-rich conditions. 



position of defect 


(eV) 




4.87 




4.72 




4.73 



FIG. 1: (Color online) (a) The structure of the defect free ZnONW 
with marks of several positions for doping, (b) The ZnONW of two 
periodic units with Ga, N and P dopants at bl, bl and 62 positions, 
respectively, represented by ZnONW:(Ga''^+N''^,P''^). The atoms of 
Zn, O, N, P and Ga are colored in gray, red, blue, green and purple, 
respectively, throughout the context. 



and the parameter C of the defect can be obtained by calcu- 
lating E^(D'^,l) with an axial length I of different unit cells, 
and then linearly fitting^" E^{D'',l) to the inverse of I. Thus, 
the error of the defect formation energy AEf{D'^,l) can be 
evaluated by 

AEf{D\l)^El{D'^)-Ef{D\l). (7) 

The supercells of ZnO bulk, ZnO film and ZnONW con- 
tain 128, 144, 96 atoms in total, respectively. The ZnO film 
is modeled perpendicular to the axis along the [1120] direc- 
tion, containing 9 central-symmetric layers, 13.3 A thick. The 
central three layers are fixed to the bulk positions while the 
top and bottom three layers are allowed to relax without con- 
straints. The supercell of ZnONW contains two periodic units 
along the [0001] direction with diameter of 9.9 A, and the ax- 
ial length of the supercell is 10.84 A. The vacuum distances 
between the surfaces of neighboring ZnO film and ZnONW 
are kept 9.7 and 15.1 A, respectively. The size of the supercell 
of ZnONW is 25.0 x 25.0 x 10.84 A. The lengths of these three 
supercells of ZnO bulk, ZnO film and ZnONW are marked by 
2x X 2y X 2z, v x 2y x 2z and i; x i; x 2z, respectively. For 
each case, x, y and z represent the unit lengths of the axial di- 
rections, and V represents the given length of supercell in the 
direction with vacuum. By default, these supercells are used 
to calculate the formation energies E^ {D'^,l=2z) directly by 
Eq. (2). 

To evaluate the error of the defect formation energy 
A£-^ {D'i,2z) by Eq. (7) and the linear fitting method^, the 
total energies of a a; x y x z supercell of ZnO bulk, a w x y x 2; 
supercell of ZnO film, v x v x z and v x v xiz supercells of 
ZnONW are also calculated. The influence of the energy error 
AEf{D'^,2z) will be analyzed in the context. 



III. Negative Charged States of N-doped ZnO Nanowires 

The structure of the defect free ZnONW is shown in Fig. 
1(a), where several doping positions are marked. To investi- 
gate the properties of the defects in the supercell of ZnONW 
with two periodic units, for every position only one of the 
two units is substitutionally doped in order to avoid two 



TABLE IL The formation energies E-^ of neutral charged 
ZnONW;2N under O- and N-rich conditions. 

positions of defects E^ieV) 





9.37 


2^61=2 


9.39 




9.52 




9.61 


2jyc2c3 


9.53 


2^c2c7 


9.55 



same dopants occupying the same position. As an exam- 
ple, ZnONW with Ga, N and P dopants at 61, 61 and 62 
positions, respectively, is shown in Fig. 1(b), marked as 
ZnONW:(Ga"H-N",P''2), where (Ga^H-N") represents that a 
pair of Ga and N dopants at 61 are bonded with each other Ga 
substitutes Zn atom, N and P substitute O atoms, respectively. 

To compare the stability of the structures with different de- 
fects, the formation energies are calculated. For the defect 
No in the ZnONW supercell (ZnONW:N), all a, 6 and c posi- 
tions are equivalent, respectively. The formation energies are 
Hsted in Table I. The difference of E^ between ZnONW:N^i 
and ZnONWiN'^i is 0.01 eV, while E^ of ZnONW:N''i and 
ZnONWiN'^i are lower than that of ZnONW:N"i by at least 
0.14 eV. For all of these cases the total magnetic moment per 
supercell after full relaxation is found to be 1.000 /is. It indi- 
cates that the surface (c) and subsurface (6) sites are energet- 
ically favorable for the Nq defect. For doping two N atoms 
in the supercell (ZnONW:2N), E^ of several selected cases 
are shown in Table II. Among these cases, ZnONW:2N''^'^^ is 
energetically the most favorable. These results are in general 
agreement with the results of a previous work^'. In the fol- 
lowing, only ZnONW:N''i and ZnONW:2N^i''i are discussed 
as typical examples with one or two doped N atoms in the 
ZnONW, if not specified. For a few interesting cases, more 
than one configurations of energetically favorable defects with 
nearly the same E-^ will be considered. 

Based on Eq. (2), E^' ofZnONW:N''i andZnONW:2N''i'=i, 
and doping one N atom in the supercells of ZnO bulk and film 
in different charged states are shown in Figs. 2(a) and (b), 
(c) and (d), respectively, where the negative range under the 
Fermi level (FL) corresponds to the valence band. As the FL 
moves from the valence band maximum (VBM) to the con- 
duction band minimum (CBM), the semiconductor changes 
from p-type to n-type. In Fig. 2, as the band-gap gives the 
upper limit of the FL, the band-gap of ZnO bulk is taken as 
the experimental value of 3.43 eV-*^, while the band-gaps of 
ZnO film and ZnONW are approximately set to 3.45 and 3.70 
eV, respectively. 
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FIG. 2: (Color online) The formation energy of defect No in 
ZnONW (a) under O- and N-rich, (b) under Zn- and N-rich con- 
ditions, in ZnO bulk and film (c) under O- and N-rich, and (d) under 
Zn- and N-rich conditions. The cases of ZnONW:N, ZnONW:2N, 
ZnO film:N and ZnO bulk:N are colored in red, blue, red and black, 
respectively. Relatively stable states above the VBM are drawn in 
solid (ZnO NW/film:N) or dashed (ZnONW:2N/ZnO bulk:N) lines, 
while the unstable states above the VBM and all states below the 
VBM are indicated in dots, throughout the context. 



For ZnO bulk, the (0/-1) charge transition of Nq defect is 
deep in the gap, in agreement with the previous observation^^. 
On the other hand, the (0/-1) charge transition of ZnO film 
happens below the VBM; more than five transitions to nega- 
tive charged states of ZnONW appear within the valence band. 
It is reported both experimentally and theoretically that the 
band of surface O 2p dangling bond of ZnO sits below the 
VBM.^^'^'' This may be true for low dimensional ZnO but not 
for all metal-oxide semiconductors^. So, the charge transi- 
tion points below VBM in Fig. 2 correspond to the unoccu- 
pied surface states in the ultrathin ZnO film and ZnONW. As 
a result, in contrary to bulk, Nq in ZnO film and ZnONW 
have many negative charged states in the band-gap, especially 
under n-type condition with the FL close to the CBM. These 
negative charged states are easier to get negative E-^ than ZnO 
bulk. 

However, not all of these negative charged states are stable. 
For ZnONW and n-type ZnO film, as the defect carries more 
negative charge, the variation of for adding one more neg- 



ative charge becomes smaller rapidly. When the energy dif- 
ference between different charged states is smaller than the 
energy lowering for producing a Zn vacancy (Vzn) with the 
same negative charge, it will induce spontaneous disintegra- 
tion from higher negative states into lower negative states and 
Vzn defects^. The most stable state of Vzn is -2 charged 
(Yzn) ill spi*^^ °f O-rich, Zn-rich, p-type and n-type condi- 
tions. For ZnO film, Ef {Vz;,) is -0.76, -7.66, 2.81 and -4.09 
eV for O-rich p-type, O-rich n-type, Zn-rich p-type and Zn- 
rich n-type extremes, respectively. For ZnONW, these values 
are -6.48, -13.86, -3.13, and -10.51 eV, respectively. An ap- 
proximate criterion to filter the unstable charged states of Nq 
can be written as 



Ef{N^o) + -Ef{Vl-) 



< Ef{N^o^) (q < 0). (8) 



If the above inequality is satisfied, the {q — 1) charged state 
will be unstable; otherwise the charged states are relatively 
stable without self-disintegration, as shown in Fig. 2. 

For ZnONWiN" and ZnONW:2N"'=\ as manifested in 
Figs. 2(a) and (b), the relatively stable states free from the 
spontaneous disintegration include the neutral and -1 charged 
states under the O-rich condition, and 0^-2 charged states un- 
der Zn-rich condition, among which the ground state in energy 
is - 1 charged under O-rich condition and -2 charged under Zn- 
rich condition, respectively. In general, E^ under the O-rich 
condition is larger than that under the Zn-rich condition and, 
it is easier to get a ground state with a higher negative charge 
under the Zn-rich condition than that under the O-rich condi- 
tion. It is evident that ZnONW:2N''i'=i has a larger E^' than 
ZnONW:N". 

The negative charged states with negative E^ play roles as 
acceptors to compensate electrons, resulting in high compen- 
sation rate and low photocurrent. The equilibrium concentra- 
tion of a defect D at temperature T is given by25i2^ 



CiD,q) a exp[-Ef{D,q)/kBT]. 



(9) 



Obviously, a larger positive E^ leads to a lower defect con- 
centration, and for a negative E-^ , the defect can be viewed 
as an active compensation center. So, a small positive E^ 
of the ground state is prefeiTed, which leads to a high de- 
fect concentration and a low recombination rate. To meet 
this requirement, among the ground states of ZnONW:N''^ and 
ZnONW:2N''^'=\ E^ under the O- and N-rich environments of 
ZnONW:(N''i)i- at p-type extreme and ZnONW:(2N''i'=i)i- 
at n-type extreme are favorable, whose E-^ are 0.67 and 1.73 
eV, respectively. 

The defect energy errors of the charged states are eval- 
uated by Eqs. (5) and (7) with the linear fitting method. 
AE^ [(Nq )i^,2z] is found to be -0.43 eV, and the correspond- 
ing parameter C is 4.94 eV A. It means that after the energy 
coii-ection, in Fig. 2(a), E^' of ZnONW:(N''i)i- under O- and 
N-rich environments at p-type extreme drops from 0.67 eV 
to 0.24 eV. Moreover, AE^ [(N'^^^f- ,2z] is -1.48 eV, with the 
parameter C of 15.50 eV-A; AEf[(Wzn?^,2z] is -1.41 eV, 
with C of 15.27 eV-A. Substituting these three values of AE^ 
into (8), this inequality is still satisfied. So, under the O- and 
N-rich environments, the -1 charged state of ZnONWiN**^ is 



still the ground state, and at p-type extreme of 0.24 eV 
will also lead to a high defect concentration and a low recom- 
bination rate. It can be seen that as the negative charge in- 
creases from -1 to -2, AE^ of ZnONW:N''^ is enlarged from 
-0.43 eV to -1.48 eV. For the -3 charged state of ZnONW:N''i, 
E^ is expected to be corrected downward by more than 1 .48 
eV, and it may turn from the unstable state to stable state un- 
der the Zn- and N-rich environments in Fig. 2(b). However, 
the -3 charged state is less important and has no influence on 
the water splitting. Similarly, after an energy correction, al- 
though E^' of ZnONW:(2N''i'=i)i- is lowered by a fraction 
of one electi-on volt, ZnONW:(2N''i^i)i- under O- and N- 
rich environments at n-type extreme is still the only state of 
ZnONW:2N''^'=^ favorable for the water splitting. 

As a comparison, the charged states of ultrathin ZnO film 
are also discriminated by the same criterion in inequality (8), 
as shown in Figs. 2(c) and (d). Except for the narrow range 
near the n-type extreme under the O-rich environment, the 
ground state of No defect in ZnO film is highly negative 
charged with a negative E^, which is easy to lead to a high 
recombination rate. If E-^ can be lifted from negative to pos- 
itive, it would suppress the recombination rate and raise the 
photocurrent. Normally, the passivated codoping could lead 
to a E^ larger than the single doping, as evidenced by the ex- 
perimental observation that (Ga-nN) codoped ZnO films have a 
significantly enhanced photocurrent and a suppressed recom- 
bination rate than the solely N doped ZnO film^. The present 
analysis is effective to explain the experimental phenomena 
for both ZnONW and ZnO film. For the -1 charged state, 
AE^ from Eq. (7) is -0.50 and 0.21 eV for ZnO film and 
bulk, respectively. However, such energy errors will not alter 
the conclusions for ZnO film and bulk. 

Fig. 3 shows the DOS and optical absorption spectra of 
ZnONWiN" and ZnONW:2N''i'=i. For bofli cases, the par- 
tial DOS (PDOS) of N in Figs. 3(a) and (b) reveals that the 
valence electrons of N dopants accumulate at the top of the 
valence band, resulting in a group of defect bands above the 
VBM of pure ZnONW. Compared with the pure ZnONW, the 
band-gaps of ZnONW:N" and ZnONW:2N^i'=i are narrowed 
by 0.16 and 0.24 eV, respectively. 

Figs. 3(c) and (d) present the optical absorption spectra 
of ZnONWiN" and ZnONW:2N''i^i with different charged 
states, respectively, in comparison with the pure ZnONW. 
It is noticeable that even for the negative charged states in 
Figs. 3(c) and (d), the defect bands in the band-gap are not 
so close to the CBM as the gap is smaller than 0.5 eV. In 
fact, the spectra of the negative charged states in the range 
of 0.0-0.5 eV are contributed by the electron transition from 
the main valence bands to the defect bands, or within the de- 
fect bands between the VBM and CBM. So are the spectra of 
neutral ZnONW:(N''i)0 from 0.00 to 1.30 eV, and the spec- 
tra of neutral ZnONW:(2N''i'=i)0 from 0.00 to 1.22 eV. For 
the application in the PEC water splitting, only the electron 
transition from the valence bands into the conduction bands 
is meaningful and effective. Thus, such spectra in the range 
close to 0.0 eV are ineffective, and will not be discussed. Fo- 
cusing on the remaining part of the spectra, the spectra of - 
1 charged ZnONW:N^\ -1 and -2 charged ZnONW:2N''i=i 
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FIG. 3: (Color online) The DOS of the neutral state of No doped 
ZnONW (red solid line): (a) ZnONW:N" and (b) ZnONW:2N"=\ 
in comparison to the pure ZnONW (black dash line). The partial 
DOS of N (blue filled area) dopants are also included. The zero point 
of energy is set to the VBM of pure ZnONW. The Fermi level with N 
dopants are indicated in red short dash line. The band-gap values are 
marked. The optical absorption spectra of No '■ (c) ZnONWiN*'^ and 
(d) ZnONW:2N"''\ The spectrum of pure ZnONW is indicated by 
the olive dash-dot line, the spectra of 0, -1 and -2 charged states are 
plotted with black dashed line, red solid line and blue dots, respec- 
tively. 



have red-shifts more remarkable than those of either their neu- 
tral charged states or the pure ZnONW. A common ground of 
these low negative charged states is that their total charges are 



no larger than that of the corresponding defect-free nanowires. 

Moreover, by incorporating the results from Figs. 2(a) 
and (b), we find that ZnONW:(N''^)^^ at p-type extreme and 
ZnONW:(2N''i^i)i" at n-type extreme under the O- and N- 
rich environments are the two most probable states of the ul- 
trathin ZnONW with Nq defects to perform a nice PEC re- 
sponse. These two states are ground states with proper for- 
mation energies on one hand, and have the optical absorption 
spectra with significant red-shifts on the other hand. If the 
spectra of neutral charged ZnONWs remain in the range of 
ultraviolet light, only these negative charged states are pos- 
sible to red-shift the spectra into the visible range, leaving a 
series of absorption peaks from the ultraviolet to visible range. 

It is intriguing to note that for N doped n-type ZnONW 
for PEC water splitting, the optical absorption spectrum was 
experimentally observed to leave a long decaying tail from 
the ultraviolet to visible range, where the tail absorbance at 
500 nm is approximately the half of the main absorption peak 
at 400 nm". In experiments it was identified that the domi- 
nating defect in ZnONW is Nq, with the N concentration of 
3.7%". Comparing this experiment to our theoretical results 
of ZnONW:N" and ZnONW:2N''i'=\ we deduce that the ex- 
perimentally observed spectrum with a tail stretching into the 
visible range can be attributed to the negative charged ground 
states of No- 



IV. Uncompensated Codoping of (Ga, N, P) in ZnO Nanowires 

To consider the case of Pq defect in ZnONW (ZnONW:P), 
as the position c is energetically favorable for Pq, 
ZnONWiP"^"^ is chosen as an example to investigate the gen- 
eral properties of ZnONW:P. Fig. 4(a) shows that the defect 
bands of P situate deeper than N in the band-gap of pure 
ZnONW, giving rise to a sharp peak of DOS, and the band- 
gap is narrowed from 1 .46 eV of the pure ZnONW to 1.16 eV. 
Aiming at narrowing further the band-gap, it is conceived that 
N and P are codoped into ZnONWs. Combining the charac- 
teristics of N and P dopants, the band-gap may be narrowed 
remarkably by P while the DOS of the defect bands is main- 
tained to spread from the VBM as continuously as possible by 
either N or P dopants. Accordingly, by incorporating the ad- 
vantages of passivated codoping by keeping the crystallinity, 
controlling the formation energy, and suppressing the recom- 
bination rate^iS, we try to codope one (Ga-nN) and one (Ga-nP) 
in pair per supercell into ZnONWs [ZnONW: (GaH-N,GaH-P)]. 

To examine the above consideration, as an example, the 
DOS, formation energy and optical absorption spectrum of 
ZnONW:(Ga"H-N",Ga'"'H-P'"*) are calculated, as presented in 
Figs. 4(b)-(e). From Fig. 4(b), one may see that the band- 
gap of ZnONW: (Ga^H-N" ,Ga^4+pf,4-) narrowed to 1 .08 eV, 
which is smaller than that of ZnONW:P'^'^. Meanwhile, the 
gap between the deepest P defect bands and the lower defect 
bands of N and P is kept quite narrow. 

could be manipulated by the chemical potential of 
P (/ip) between the values of P-rich and P-poor^^, where 
it is known that the most P-rich potential of is calcu- 
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FIG. 4: (Color online) The DOS of the neutral state for 
(a) ZnONW:?'^^ and (b) ZnONW:(Ga"-^N",Ga'"'-HP*"'), includ- 
ing the PDOS of N (blue) and P (green), which are com- 
pared with that of pure ZnONW. The formation energy of 
ZnONW: (Ga"-HN",Ga"-^P'"') including P-rich (black) and P-poor 
(red) cases under (c) 0-, N- and Ga-rich, and (d) Zn-, N- and Ga-rich 
conditions, respectively, (e) The optical absorption spectra of pure 
ZnONW and the charged states of ZnONW:(Ga"-HN",Ga'"'-^P''''). 



lated from P4O1G , and we choose P4 to obtain a P-poor /ip approximately. Figs. 4(c) and (d) show that under 



TABLE III: The formation energies of neutral charged 

ZnONW:(Ga,N,P) under 0-, Ga-, P-, and N-rich conditions. 

positions of defects 'W\eVY 

{Ga"' +iV''\P"^) 15.78 

(Ga" + iV",P''2) 15.75 

(Ga" + Af^.P*"*) 15.80 

(Ga" + Ar^.P*"^) 15.78 

(Go" +P",iV''^) 16.2 

{Ga''^ + P''^ , N"'^) 16.2 

(Ga" + P",iV'"') 16.2 

(Ga"2 + P''^iV") 20.1 

(Ga",Ar",P''2) 

(Ga" +iV" + P''i) 16.1 

(Ga" + iV",P°^) 16.3 

(Ga"' + AT" + P"2) 16.2 

{Ga"'^ + N'''^,P''^) 16.3 

(Ga"'^ + iV''^ P''^) 16.4 



the O-rich condition none of the negative charged states of 
ZnONW:(Ga"+N^\Ga'"*+P''4) is stable, and the ground state 
is neutral charged [ZnONW:(Ga"+N",Ga'"*+P''4)°] with a 
positive , whereas under the Zn-rich condition the ground 
state is -2 charged [ZnONW:(Ga"+N",Ga''4+P''4)2-] with 
a negative at n-type extreme. However, in Fig. 4(e), 
the optical absorption spectra of -1 and -2 charged states 
have no effective red-shifts from the neutral charged state 
that can narrow the gap between the defect bands and 
the CBM. This is because either -1 or -2 charged state 
of ZnONW:(Ga''i+N''i,Ga''4+P''4) has a total charge larger 
than that of the defect-free nanowire. Therefore, no matter 
whether the formation energy is corrected by Eq. (7) or not, 
ZnONW:(GaH-N,GaH-P) is not so useful for the water splitting. 

To improve the performance, it turns out to take advan- 
tage of the low negative charged ground state with a total 
charge equal to the defect-free nanowire, as is in the case 
of ZnONW:N''^. A strategy can be proposed to promote 
the PEC responses of ZnONWs by means of the uncom- 
pensated codoping of one (Ga-nN) in pair and one single 
P dopant per supercell [ZnONW:(GaH-N,P)]. E^ for several 
ZnONW:(Ga,N,P) cases are listed in Table III. It is seen that 
the defects including a pair of (Ga^H-N^i) and one single P at 
a2 or 62 or 63 or 64 position are relatively the most favorable 
configurations in energy. The difference of E^ between these 
four cases are so small that it is no larger than 0.05 eV. Com- 
paring the cases of (Ga"H-N^\P) and (Ga''iH-P''\N), one finds 
that the pair of (Ga^^H-N''^) is energetically more favorable 
than (Ga^i-nP"). In presence of the (Ga^H-N^i) pair, the P 
dopants at c positions are unstable, because E-^ in these cases 
is usually larger than 200 eV. The three typical examples such 
as ZnONW:(Ga"H-N''i,P''2), ZnONW:(Ga''iH-N''i,P'^2) ^^^j 

ZnONW:(Ga"H-N",P'"*) are energetically favorable, whose 
total magnetic moments per supercell after full relaxation are 
found to be 1.093, 1.158 and 1.146 /is, respectively. 

As shown in Figs. 5, all the band-gaps of these 
three cases are narrower than those of ZnONW:N^^, 
ZnONW:2N"'=i and the pure ZnONW, but are close to 
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FIG. 5: (Color online) The DOS of the neutral state for (a) 
ZnONW:(Ga"+N",P''^), (b) ZnONW:(Ga"+N",P''2) and (c) 
ZnONW:(Ga"+N",P''*), including the PDOS of N and P which are 
compared with that of pure ZnONW. 



the band-gap value of ZnONW:(Ga"H-N",Ga'"*H-P''4). por 
ZnONW:(Ga^iH-N",P'^2) ZnONW:(Ga^iH-N^i,P''4), there 
are also defect bands spreading between the VBM and the 
deepest defect peak of DOS. 

In Figs. 6(a) and 6(b), the ground states of 
ZnONW:(Ga^iH-N''\P''2) under the O-rich and Zn-rich en- 
vironments are -1 [ZnONW:(Ga''iH-N",P''2)i-] and -2 
[ZnONW:(Ga"H-N^i,P''2)2-] charged, respectively. E^' of 
ZnONW:(Ga^iH-N''i,P''2)i- is -0.54 eV at n-type extreme 
under the 0-, N-, Ga-rich and P-poor environments [Fig. 
6(a)]; Ef of ZnONW:(Ga^iH-N*i,P''2)2- is -0.94 eV at n- 
type extreme under the Zn-, N-, Ga- and P-rich environments 
[Fig. 6(b)]. For these two states, evaluating the forma- 
tion energy error from Eq. (7), AEf[(Ga'''^+N''\P''^)^- ,2z] 
is -0.46 eV with the parameter C of 4.99 eV A; and 
A£;^[(Ga''iH-N^i,P''2)2-,22] is -1.53 eV with C of 15.88 
eV-A. Thus, after counting the correction on energy, E^ 
changes from -0.54 eV to -1.00 eV, and from -0.94 eV 
to -2.47 eV. No matter with or without the energy correc- 
tion, these two cases are kept to be in the ground states 
and can easily produce a small positive E^ in the ground 
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state by manipulating the n-type extreme and the chemi- 
cal potential /ip. They are energetically preferable for a 
high defect concentration and a low recombination rate. 
of ZnONW:(Ga"H-N",P'^2) and ZnONW:(Ga"H-N''i,P''4) are 
similar to that of ZnONW:(Ga"H-N''\P''^) because of a very 
small Ef difference among them. 

Fig. 6(c) manifests that for ZnONW:(Ga^iH-N^\P'"*)i-, 
ZnONW:(Ga"H-N",P''2)i- and ZnONW:(Ga"H-N''\P'^2-)i-^ 

the leftmost absorption peak with an intensity larger 
than 1.0 Hes at 0.72, 0.75 and 1.00 eV, respectively. 
Among them, the spectra of ZnONW:(Ga"H-N",P''4)i- and 
ZnONW:(Ga"H-N",P''2)i- have further red-shift from those 
of ZnONW:(N")i- and ZnONW:(2N^i'=i)i^. The peaks 
around 1.62 eV correspond to the band-gap from the main 
VBM to the main CBM, in spite of the defect bands. In con- 
trast, ZnONW:(Ga''^H-N''\P''2)2" gives no efficient band-gap 
narrowing, so are the -2 charged states of the other two cases. 

Due to the severe underestimation of band-gaps by the DFT, 
all the results above calculated within the GGA cannot quan- 
titatively reflect the realistic situation in PEC responses. For 
this purpose, the GGAh-U and extrapolation methods that are 
introduced in Section II are employed here. When U correc- 
tion of {U-J)=5.Q eV for the d electrons is added to the GGA 
calculation, the band-gaps and spectra are partially corrected. 



as shown in Fig. 7(a). For ZnONW:(Ga''iH-N",P'"')i-, 
ZnONW:(Ga^iH-N^\P''2)i-, ZnONW:(Ga''iH-N''i,P°2)i-, 
ZnONW:(N''i)i- and ZnONW:(2N''i'=i)i-, the leftmost 
absorption peak with an intensity larger than 1.0 lies at 0.93, 
0.96, 1.22, 1.29 and 1.42 eV, respectively. 

From the preceding results, we unclose that under 
the 0-, N-, Ga-rich, P-poor and n-type environments, 
the uncompensated (Ga, N, P) codoping cases such 
as ZnONW:(Ga"H-N''i,P''4)i-, ZnONW:(Ga"H-N''i,P''2)i- 
and ZnONW:(Ga"H-N",P''2)i- are energetically favorable, 
whose total charge is equal to that of the defect-free nanowire. 
In general, their optical absorption spectra have larger red- 
shifts and higher absorption peaks than the cases of exclusive 
N doping. Meanwhile, the (Ga, N, P) codoping cases have 
proper formation energy that can lead to a high defect concen- 
tration and a low recombination rate. Both factors are much 
favored in enhancing the photocurrent in the PEC water split- 
ting. 

To show the improvement of the uncompensated codoping 
of (Ga, N, P) compared with the exclusive N dopoing in the 
ZnONWs, in a way approaching to the experiment, an approx- 
imate estimation is performed by utilizing the extrapolation 
equation (1). The calculated results by GGA and GGAh-U 
with a ([/ — J) of 5.0 eV for d electrons can be extrapolated to 
the experimental ones. ZnONW:(N''^)^^ is taken as the stan- 
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dard defect of Nq in Eq. (1), and E'^^p is about 3. 1 1 eV, which 
corresponds to the red-shifted peak of the N doped ZnONW 
at about 400 nm in the UV-vis spectra in the experiment". 
Consequently, the resuhed Eg of ZnONW:(Ga"+N''\P''4)i-, 
ZnONW:(Ga"+N''i,P''2)i- and ZnONW:(Ga"+N''i,P"2)i- 
is 1.96, 1.99 and 2.30 eV, coiTesponding to 633, 623 and 540 
nm, respectively. 

Besides the extrapolation method, to give another direct and 
intuitive illustration, U correction can be added for both d 
and p electrons. When {U — J) is assumed as 18.0 eV for 
d electrons and 15.0 eV for p electrons, respectively, the long- 
wavelength end of the optical absorption spectra of the exclu- 
sive N doping is manipulated to stop just at the range of violet 
and blue light, which is close to the experiment". Such a set- 
ting of GGAh-U focuses on directly correcting the band-gap 
to approach the experimental value, while the expense is that, 
the inner band structures are disturbed more greatly than the 
extrapolation method with a small U added only for d elec- 
trons. As shown in Fig. 7(b), at the long-wavelength end 
of the spectra of ZnONW:(N")i- and ZnONW:(2N"'=i)i-, 
the last peak with an absorption intensity larger than 0.1 
locates at 463 and 438 nm, respectively. In contrast, the 
last peak with an absorption intensity larger than 0.3 of 
ZnONW:(Ga^iH-N^\P'"')i-, ZnONW:(Ga^iH-N''i,P''2)i- and 
ZnONW:(Ga"H-N",P"2)i- locates at 675, 623 and 589 nm, 
respectively. Generally, the wavelengths of these three peaks 
are close to and consistent with those obtained from the ex- 
trapolation method. 

Revealed from the estimations either by the extrapolation 
method or by directly adding U for both d and p electrons, the 
absorption spectra of the exclusive N doping in the ZnONWs 
can enter the range of violet light but hard to reach the cen- 
ter of the visible range, which just manifests the case in the 
experimentii. Judging also from the above results by either of 



the two estimation methods, one may see that for both N dop- 
ing and uncompensated (Ga, N, P) codoping in ZnONWs, the 
absorption peaks from about 400 nm to larger than 600 nm 
including their tails can cover nearly the whole visible light 
range, which may thus lead to the raise of the photocurrent. 



V. Summary 

In summary, by means of the first-principles calculations 
within the density functional theory, the charged states and 
band-gap narrowing in different doped and codoped ZnONWs 
are investigated and compared with those of ZnO bulk and 
film. It is found that the charged states of nitrogen related de- 
fect, particularly the low negative charged ground state with a 
total charge no larger than the defect-free nanowire, play a key 
role for the PEC response, which leads to a long tail in the op- 
tical absorption spectrum. A design of uncompensated codop- 
ing of N, P and Ga in ZnONWs is proposed, which could give 
an evident enhancement of the PEC response. These observa- 
tions will in general be useful for improving the PEC response 
particularly for the water splitting of wide-band-gap semicon- 
ductors. 
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